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ABSTRACT 

Absolute  total  and  thermal  neutron  flux  of  the  U,  So  Naval  Post¬ 
graduate  School's  AGN-201  reactor  was  determined  by  neutron  activation 
of  thin  gold  foils .  Foil  activities  were  measured  with  a  gamma-ray 
scintillation  spectrometer,  using  methods  designed  to  minimize  the 
effect  of  changes  in  spectrometer  gain.  Flux  values  were  calculated 
for  nominal  power  levels  of  0.1  watt  and  1,  10,  100,  and  750  watts. 
Methods  and  results  are  compared  with  those  of  previous  investigations . 
The  flux  level  was  found  to  be  a  linear  function  of  power  within  this 
range;  total  and  thermal  average  fluxes  were  determined  to  be  respectively 
6.64  X  10^  and  5.41  x  10^  neutrons  per  square  centimeter  per  second  per 


watt. 
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1.  Introduction. 

Activation  and  decay  measurement  of  gold  foils  has  become  a  standard 
technique  for  determining  reactor  neutron  flux.  For  the  U.  S.  Naval  Post¬ 
graduate  School's  AGN-2 01  reactor,  this  method  has  been  used  by  Kelly  and 
Clements  (10)  to  determine  the  absolute  thermal  flux  at  0.1  watt,  by  Ferguson 
and  Harvey  (6)  to  determine  relative  flux  and  flux  distribution  at  several 
power  levels ,  and  by  Copeland  and  Reasonover  (4)  to  determine  the  flux 
perturbation  caused  by  the  presence  of  the  foil.  In  the  present  investiga¬ 
tion,  the  absolute  flux  was  measured  at  several  power  levels,  from  0.1 
watt  to  750  watts. 

For  measuring  the  absolute  disintegration  rate  of  the  irradiated  foil, 
several  methods  have  been  developed.  Those  which  depend  on  beta-count¬ 
ing  require  the  least  special  equipment,  but  they  involve  the  corrections 
and  difficulties  associated  with  absolute  beta  measurements.  It  is  possible 
to  compare  the  activation  induced  by  the  reactor  with  that  from  a  standard 
neutron  source;  this  technique  is  also  easy  to  apply,  but  the  results  are  only 
as  good  as  the  accuracy  to  which  the  neutron  density  and  energy  distribution 
of  the  source  are  known.  For  speed  and  overall  accuracy,  gamma-ray 
spectrometry  presents  several  advantages,  and  spectrometry  was  the  tech¬ 
nique  used  in  the  present  investigation.  The  characteristics  of  scintillation 
spectrometers  require  precautions  to  be  taken  against  drift  of  the  photo¬ 
multiplier  tube  and  of  the  counting  circuits  themselves;  this  drift  occurs 
from  several  causes,  and  it  is  not  always  apparent  from  the  results  of  a 
short  series  of  counts.  In  order  to  compensate  for  the  effects  of  drift,  a 
somewhat  novel  procedure  was  used  for  determining  the  parameters  of  the 
principal  gold  photopeak ,  that  at  0.411  MEV. 
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2.  Experimental  procedure. 

Thin  circular  gold  foils,  with  nominal  dimensions  of  0,5"  x  0,0005" 
and  a  mass  of  about  30  mg,  were  weighed  to  the  nearest  0,5  mg  and  were 
rapidly  inserted  into  the  glory-hole  of  the  reactor  after  it  had  been  stabilized 
at  the  desired  power.  The  foils  were  located  at  the  center  of  the  core  to 
an  estimated  accuracy  of  +  1  mm,  and  insertion  and  removal  times  were 
controlled  to  +  1  second.  For  each  run,  the  sample  holder  also  contained 
an  indium  foil  monitor,  which  was  located  in  the  graphite  reflector  9" 
from  core  center.  This  location  was  chosen  to  minimize  flux  depression 
from  the  indium,  while  still  exposing  it  to  a  significant  neutron  population. 

In  addition,  the  characteristics  of  the  reactor  are  such  that  the  flux 
distribution  in  the  reflector  is  relatively  constant  compared  to  that  at  the 
edges  of  the  core;  thus  any  position  error  would  have  minimum  effect. 

The  results  of  Ferguson  and  Harvey  showed  that  at  the  location  of  the 
monitor  foil  the  epithermal  flux  is  negligible ,  so  that  no  correction  for 
fast  flux  was  required;  this  conclusion  may  not  be  valid  for  power  levels 
above  100  watts. 

Foils  were  irradiated  at  0,1,  1,  10,  100,  and  750  watts.  Although 
the  reactor  can  be  brought  to  1000  watts  for  a  short  time,  it  was  difficult 
to  maintain  accurate  control  as  the  sample  holder  was  inserted  and  to 
reproduce  conditions  exactly  for  duplicate  runs.  Insertion  of  the  sample 
holder  causes  an  unavoidable  change  in  reactivity,  which  requires  adjust¬ 
ment  of  the  control  rods,  and  at  high  power  levels  the  uncertainty  in  the 
neutron  flux  to  which  the  foil  was  exposed  and  in  the  timing  becomes 
greater.  To  minimize  these  errors  would  have  required  an  exposure  time 
that  would  have  produced,  at  1000  watts,  an  unacceptably  high  level  of 
activity,,  because  of  the  length  of  time  the  samples  would  have  had  to  de¬ 
cay  before  counting.  For  these  reasons,  750  watts  was  the  maximum 
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level  at  which  measurements  were  attempted. 

Duplicate  runs  were  made  on  all  samples.  Irradiation  time  s  were 
the  same  for  both  runs  at  each  power  level ,  except  that  during  the  first 
activation  at  0 .  Iw  a  line  voltage  transient  caused  an  undesired  scram  which 
interrupted  the  run.  The  irradiation  time  given  is  corrected  for  sample 
decay  during  the  interruption.  The  exposure  times  for  the  various  runs 
were: 

0.1  watt  —  477.7  minutes  (corrected)  and  534  minutes 

1  watt  —  60  minutes 

10  watts  —  30  minutes 

100  watts  —  10  minutes 

750  watts  —  5  minutes 

Several  factors  affected  the  choice  of  irradiation  time.  A  minimum 
of  ten  minutes  was  considered  desirable,  in  order  to  minimize  the 
relative  error  caused  by  in sertion  and  removal  times.  For  the  lower 
powers,  additional  time  was  required  to  bring  the  activity  to  a  level  at 
which  the  photopeak  maximum  would  be  at  least  the  several  thousand 
counts  per  minute  necessary  to  give  a  clear  and  sizable  photopeak  and 
to  reduce  the  uncertainties  produced  by  counting  statistics  and  by 
background.  On  the  other  hand,  to  prevent  coincidence  losses  in  the 
spectrometer  the  maximum  count  was  not  allowed  to  rise  above  about 
^20,000  cpm;  this  meant  that  the  100  watt  samples  had  to  decay  for  about 
six  days  before  being  counted,  and  the  effect  of  uncertainty  in  the  accepted 
2.7  day  half-life  of  Au-19 8  could  be  significant.  The  5 -minute  time 
chosen  for  the  750  watt  run  was  a  compromise  between  decay  and  timing 
errors . 
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Upon  removal  from  the  reactor,  the  indium  and,  when  necessary,  the 
gold  foils  were  allowed  to  decay  until  their  activities  were  at  a  suitable 
level  for  counting.  The  indium  was  counted  in  a  standard  G-M  counter; 
the  location  of  the  foil  in  the  counter  was  carefully  reproduced  for  each 
run,  but  since  only  relative  activity  was  required,  no  absorption  or  geometry 
corrections  were  applied.  Two  integrated  10 -minute  counts  were  taken  of 
each  indium  foil.  Coincidence  and  decay  corrections  were  applied  di¬ 
rectly  to  the  integrated  count  to  determine  foil  activity;  the  method  of 
determining  the  coincidence  correction  is  given  in  Appendix  I.  In  the 
higher-power  runs,  a  significant  (up  to  1000  cpm)  activity  due  to  In~114 
was  observed.  Since  this  isotope  has  a  49-day  half-life,  it  was  only 
necessary  to  allow  each  foil  to  decay  for  24  hours;  after  this  period  the 
activity  of  In-116  was  reduced  by  a  factor  of  10^,  and  the  remaining 
activity  was  from  In-114.  The  foil  was  then  counted  a  second  time,  and 
the  second  count  was  subtracted  .  from  the  first  as  "background" . 

Whether  or  not  the  activity  of  the  gold  was  high  enough  to  require 
additional  decay  time  before  counting,  the  first  foil  from  each  power 
level  was  allowed  to  remain  in  the  spectrometer  sample  mount  overnight 
before  counting ,  in  order  to  stabilize  the  photomultiplier  tube  as  far  as 
possible.  The  foil  was  placed  on  the  "3cm"  shelf  of  the  mount,  and 
the  high  voltage  of  the  spectrometer  was  set  at  1270  volts  throughout;  it 
has  been  determined  that  this  combination  gives  good  results.  The  gain 
and  bias  controls  of  the  spectrometer  were  adjusted  to  give  a  usable  peak, 
as  discussed  below.  The  0.411  MEV  photopeak  was  counted  at  least 
five  times  for  each  foil.  Since  the  only  information  desired  was  the 
parameters  of  the  photopeak  ,  no  attempt  was  made  to  determine  the 
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entire  spectrum  ^  or  to  count  more  channels  than  were  needed  to  ensure 
Inclusion  of  the  peak.  The  counts  were  corrected  for  foil  decay  before 
the  photopeak  area  was  computed;  this  approach  permitted  immediate; 
direct  comparison  of  the  results  of  duplicate  runs.  Decay-corrected  peak 
count  data  are  presented  in  Appendix  11, 

3 .  Spectrometer  stability  and  drift , 

Any  variation  in  the  overall  gain  of  a  spectrometer  will  appear  as  a 
drift  or  shift  in  the  channel  at  which  the  photopeak  maximum  appears  for 
gammas  of  a  given  energy.  The  variables  involved  have  been  discussed 
by  Altekruse  (1) ,  Covell  and  Euler  (5) ;  and  Cantarell  (2) ,  among  others. 
Briefly;  there  are  four  primary  causes  of  channel  drift:  (a)  "fatigue";  of 
the  photomultiplier;  (b)  short-term  gain  changes  in  the  photomultiplier, 
caused  by  temperature  changes,  high-voltage  fluctuations ,  mechanical 
vibrations,  etc,;  (c)  overall  gain  changes  in  the  electronic  circuitry, 
caused  by  tubevand  component  aging;  (d)  short-term  electronic  changes 
caused  by  temperature  and  voltage  transients ,  Previous  investigators 
at  this  school  apparently  concluded  that  the  spectrometer  was  "stable" 
if  the  photopeak  maximum  appeared  in  the  same  spectrometer  channel  on 
all  runs;  that  this  is  too  broad  an  assumption  is  shown  by  the  fact  that  at 
the  gain  settings  used  by  Copeland  and  Reasonover,  a  shift  of  one  full  5- 
volt  channel  corresponds  to  a  gamma  energy  change  of  0,015  MEV,  while 
in  Kelly  and  Clements'  work  a  1-channel  shift  corresponds  to  0,032  MEV 
of  energy.  Although  a  shift  or  drift  of  almost  one  channel  was  observed 
on  one  run  during  the  present  work,  most  of  the  shifts  which  occurred 
were  of  the  order  of  a  tenth  of  a  channel  width,  and  at  the  gain  setting 
used,  one  channel  corresponded  to  only  0,010  MEV,  Although  long- 
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term  channel  drift  is  not  as  serious  in  this  type  of  work  as  it  is  in  the 
analysis  of  unknown  materials,  there  is  still  an  observable  effect,  and 
rapid  shifts  caused  by  transients  can  completely  invalidate  a  run.  The 
result  of  a  shift  is  to  distort  the  apparent  shape  of  the  photopeak;  if  the 
drift  is  "down-channel",  as  most  of  it  is,  a  progressively  smaller  fraction 
of  the  "actual"  number  of  events  appears  in  the  count  taken  on  each 
successive  channel.  If  the  area  of  the  photopeak  is  then  computed  by 
the  method  given  by  Heath  (7)  (8) ,  which  involves  fitting  a  curve  to  the 
points  on  the  high  side  of  the  photopeak,  the  result  indicates  a  narrower 
photopeak,  hence  a  lower  level  of  activation,  than  is  actually  present.* 
Similarly,  an  "up-channel"  shift  gives  activation  values  which  are  too 
high.  Because  of  the  steep  sides  of  the  normal  distribution  curve  which 
contains  the  photopeak,  a  rather  slight  shift  in  gain  can  cause  a  relatively 
large  change  in  the  computed  area.  Inspection  of  the  data  of  the  previous 
investigators  shows  that  they  did  indeed  encounter  some  drift,  which  they 
attempted  to  compensate  by  averaging  the  readings  obtained  on  each 
channel  from  several  counts.  The  validity  of  this  procedure  is  questionable; 
it  will  be  discussed  below. 

Of  the  causes  of  gain  drift  previously  listed,  the  slow  aging  of 
electronic  components  was  considered  to  be  negligible  over  the  hour  or  two 
required  for  each  set  of  counts,  although  its  effect  could  easily  be  observed 


*The  effect  of  the  shift  is  to  make  the  high  side  of  the  distribution 
curve  appear  steeper  than  it  actually  is.  The  error  comes  from  the 
time  required  to  count  each  channel;  the  true  shape  of  the  photopeak 
is  unchanged,  but  it  is  moving  to  the  left  during  the  counting  interval. 
This  apparent  steepening  of  the  curve  does  not  Involve  an  increase  in 
the  resolution  of  the  system,  which  is  about  11  %  at  0.411MEV. 
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over  a  period  of  months o  Temperature,  too,  stayed  fairly  constant  for 
any  one  counting  period,  and  its  effect  on  the  electronic  circuits  was 
minimized  in  any  case  by  leaving  the  spectrometer  on,  except  for  necessary 
repairs,  throughout  the  period  of  this  worko  The  situation  with  regard  to 
electronic  and  photomultiplier  transients  was  not  so  simple.  That  these 
transients  did  occur  was  not  doubted;  sharp  changes  in  line  voltage, 
for  instance ,  if  strong  enough  to  scram  the  reactor ,  would  certainly  affect 
the  spectrometer.  In  one  case,  an  early  difficulty  with  anomalous  counts 
was  resolved  when  it  was  noticed  that  the  questionable  counts  were 
those  taken  just  before  and  after  each  hour.  The  trouble  was  ascribed  to 
the  school's  automatic  clock- setting  signal,  whose  3600  cycle  frequency 
could  quite  easily  feed  into  the  instrument,  despite  power  supply 
regulation.  Subsequent  counts  taken  near  the  end  of  an  hour  were  checked 
very  carefully  before  being  accepted. 

The  most  important  cause  of  channel  drift  is  fatigue  of  the  photo¬ 
multiplier  tube,  Cantarell  has  shown  that  fatigue  is  caused  by  polariza¬ 
tion  of  the  dynodes  after  electron  bombardment,  which  produces  an 
"insulating"  effect.  The  amount  of  fatigue  is  a  function  of  temperature 
and  high  voltage,  but  more  directly  of  count  rate  and  gamma  energy , 

A  tube  subjected  to  a  given  rate  of  scintillation  events  changes  its  gain 
over  a  period  of  hours;  this  gain  change  may  be  as  high  as  20%.  The 
rate  of  gain  change  is^ logarithmic;  in  the  present  work,  the  effect  of 
drift  was  minimized  by  leaving  the  first  sample  of  each  duplicate  pair 
in  the  scintillator  mount  overnight.  By  the  next  day  the  tube  was  on 
the  asymptotic  portion  of  its  fatigue  curve,  and  the  effect  of  the  slight 
remaining  drift  was  reduced  by  counting  across  the  photopeak  as  fast 
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as  possible  for  each  run.*  In  view  of  the  precautions  observed,  it  is 
believed  that  the  only  significant  distortion  of  the  shape  of  recorded 
photopeaks  was  due  to  the  fortuitous  combination  of  the  "vertical" 
random  errors  of  counting  statistics  and  "horizontal"  random  errors  from 
the  gain  shifts  caused  by  unpredictable  and  uncorrectible  transients. 

4 .  Photopeak  parameter  computations . 

Even  if  the  spectrometer  were  perfectly  stable,  the  count  recorded 
on  each  channel  would  be  subject  to  a  statistical  probable  error  of  the 
square  root  of  the  count.  As  has  been  mentioned,  previous  investigators 
have  averaged  successive  counts  on  each  channel,  as  one  wDuld  do  for 
total  counts  obtained  with  a  G-M  tube .  Only  a  small  change  in  system 
gain,  however,  will  change  count  rates  by  several  probable  errors. 
Transient-induced  gain  shifts,  rather  than  statistical  variations, 
were  in  fact  responsible  for  a  majority  of  the  differences  between 
different  counts  of  the  same  sample,  as  can  be  seen  from  examination  of 
the  data  in  Appendix  II.  When  four  points  are  taken  on  each  of  two 
readings,  the  statistical  probability  of  all  four  shifting  in  the  same 
direction  is  one  in  eight.  As  it  happens,  in  over  half  of  the  cases 
observed  all  four  channels  shifted  together,  giving  a  strong  indication 
that  statistics  alone  was  not  causing  the  variation.  (In  evaluating  these 
data,  one  mustk'eep  in  mind  that  the  first  channel  used  is,  in  most 
cases,  slightly  below  the  photopeak  maximum,  whereas  the  last  three 
are  all  above  it .  On  a  shift  of  the  maximum  to  the  left ,  the  count 


*  The  automatic  readout  feature  of  the  spectrometer  was  not  used, 
because  recording  a  count  manually  is  nearly  twice  as  fast. 
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observed  in  the  first  channel  will  increase^  while  the  other  three  will 
decrease  o)  If  the  peak  shifts  themselves  were  completely  random ,  and 
if  they  took  place  at  a  random  rate,  simply  taking  the  mean  of  enough 
different  counts  would  compensate  for  shifts  as  well  as  for  decay 
statistics.  This  simple  approach  was  rejected  for  two  reasons:  (1)  because 
of  tube  fatigue,  there  was  a  net  overall  drift  to  the  left;  (2)  most  of  the 
shifts  not  attributable  to  tube  fatigue,  although  they  might  drive  the 
photopeak  maximum  in  either  direction,  occurred  at  a  rate  correspond¬ 
ing  to  five  or  ten  complete  counts  over  the  photopeak.  To  average  these 
out  would  have  required  the  recording  of  thirty  or  more  counts  of  each 
sample,  and  since  each  complete  count  required  about  ten  minutes,  this 
approach  was  uneconomical . 

For  these  reasons,  it  was  decided  to  treat  each  count  over  the 
photopeak  as  a  distinct  event,  compute  the  areas  obtained  individually, 
and  average  them  at  the  end.  By  so  doing,  the  effect  of  slow  drift  was 
made  negligible,  since  it  was  well  within  statistical  variations  during 
the  four  or  five  minutes  required  to  count  four  channels.  A  very  fast 
shift,  caused  by  a  rapid  transient,  would  give  a  curve  whose  computed 
maximum  or  area  was  so  different  from  the  remainder  that  it  could  be 
identified  and  discarded. 

The  standard  method  of  computing  the  area  under  the  photopeak  is 
due  to  Heath.  A  normal  distribution  of  events  about  the  0,411  MEV 
maximum  is  assumed;  when  this  curve  is  plotted  on  a  semi-logarithmic 
scale,  the  result  is  a  parabola.  In  evaluating  experimental  data, 
observed  counts  of  channels  at  and  above  the  photopeak  maximum  are 
used  in  order  to  avoid  distortion  Introduced  by  Compton  scattering  on  the 
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low  side  of  the  peak.  A  parabola  is  fitted  to  the  natural  logs  of 
these  counts,  and  the  parameters  of  the  associated  Gaussian  curve 
are  then  determined.  In  the  method  of  computation,  however,  the 
present  procedure  differed  somewhat  from  that  previously  used.  Kelly 
and'  Clements  took  several  counts  on  each  of  three  or  four  channels  to 
determine  their  parabola  and  fitted  the  curve  to  the  normalized  average  of 
their  results.  If  each  count  over  the  photopeak  is  considered  separately, 
taking  only  three  points  from  each  cant  does  indeed  give  a  set  of  perfect 
parabolas,  but  statistical  variation  makes  them  differ  greatly  from  each 
other.  Consequently,  in  the  present  work  the  spectrometer  gain  was 
adjusted  so  that  the  maximum  point  and  high  side  of  the  photopeak  covered 
at  least  four  channels.  A  program  was  made  up  for  the  CDC  1604  computer 
which  took  the  log  functions,  fitted  a  least-squares  parabola  to  them,  then 
gave  the  peak  abscissa,  peak  ordinate,  and  area  of  the  normal  curve  so 
derived.  Because  of  the  greater  relative  variation  cfthe  smaller  counts, 
it  was  necessary  to  introduce  a  weighting  factor.  The  weight  of  each 
count  was  made  proportional  to  its  square  root;  this  procedure  made  use 
of  the  greater  relative  precision  to  be  expected  from  the  higher  counts 
without  completely  swamping  the  smaller  ones.  If  no  weighting  factor  is 
applied,  the  statistical  variation  of  the  smaller  counts  causes  the  "tail 
to  wag  the  dog",  so  to  speak, and  makes  computed  areas  differ  by  an 
excessive  amount. 

The  photopeak  areas  obtained  for  each  count,  their  mean  for  each 
foil  activation,  and  the  standard  deviation  of  the  mean,  are  listed  in 
Appendix  III.  The  precision  of  results  computed  in  this  manner  is  about 
half  that  obtained  from  simple  averaging  of  counts  on  each  channel,  but 
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one  can  at  least  feel  confident  that  no  systematic  distortion  of  the 
photopeak  is  giving  consistently  high  or  low  results. 

5.  Scintillation  crystal  efficiency. 

The  sample  mount  assembly  used  had  been  carefully  constructed  to 
give  sample  distances  (for  a  thin  mounting  shelf)  of  1^  1,  3,  5^  and  10 
cm  from  the  face  of  the  scintillation  crystal «  assuming  that  the  crystal 
was  snug  against  its  can.  Upon  a  recommendation  from  Mr.  R.  L„  Heath, 
of  the  Phillips  Petroleum  Company,  X-ray  photographs  were  taken  to  check 
this  assumption,  and  it  was  found  to  be  incorrect.  There  is  a  gap, 
partially  filled  by  what  appears  to  be  a  spring  spacer,  of  9.0  mm  between 
the  0.005"  aluminum  can  and  the  face  of  the  crystal  (Figure  1) .  When 
this  distance  is  added  to  the  thickness  of  the  mounting  shelf,  a  sample 
in  the  "3  cm"  position  is  actually  4.03  cm  from  the  face  of  the  crystal. 
Efficiencies  previously  used  for  this  crystal  and  mount  had  been  taken 
from  Heath's  standard  catalogue  and  handbook  of  scintillation  spectrometry 
(7)  (8);  they  are  on  the  order  of  0.118,  for  a  3cm  distance.  The  true 
value  for  4.03  centimeters  was  computed  by  Heath  (9)  for  this  investigation; 
it  is  0.0846. 

6.  Relative  flux  from  indium  monitor  activation. 

Relative  flux  levels  derived  from  beta  counts  of  the  indium  monitor 
foils,  normalized  to  1  watt,  are  plotted  in  Figure  2.  Each  point  represents 
the  average  of  two  counts  on  each  of  two  duplicate  runs;  deviations  are 
too  small  to  plot.  As  can  be  seen,  from  0,1  to  100  watts  the  flux  is 
linear  with  power  to  the  precision  of  the  measurements.  No  explanation 
is  offered  for  the  high  values  found  at  750  watts;  the  small  deviation 
obtained  in  independent  measurements  suggests  that  this  is  a  true  value 
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Crystal  and  sample  mount 


and  not  the  result  of  random  error„  A  possible  reason  may  he  that  at. 
powers  above  100  watts  enough  epithermal  neutrons  reach  the  position  of  the 
indium  foil  to  give  a  significant  amount  of  activation  from  resonance 
capture  o 

7 .  Calculation  of  absolute  neutron  flux. 

The  number  of  events  under  the  photopeak  per  unit  time  is  related  to 
the  absolute  gamma  emission  rate  of  the  foil  by  the  expression: 


Rpt*  Et*  Fg. 


Np  =  total  (computed)  number  of  events  under  the  photopeak 
Rpt  “  peak-to-total  ratio  (0.725) 

-  crystal  detector  efficiency  (0„0846) 

Fs  =  correction  for  gamma  self-absorption  (0.997) 

F|j,  =  correction  for  internal  conversion  (0„96) 

Fg  =  correction  for  absorbing  material  in  can  (0.99) 

The  crystal  efficiency  was  provided  by  Heath  (9) ,  The  self-absorp¬ 
tion  expression  was  determined  from  the  equation: 

/at 

p_  =  - ^ - f - 

1  -  e-^t 

For  gold  foils  0.0005"  thick,  p  =  0.19,  and  t  =  0.021  g/cm^;  Fg  =  0.997. 
^ic  given  by  Raffle  (11)  as  0.96 .  F^  comes  from  the  usual  exponential 
attenuation  formula,  using  p  =  0.0287  for  aluminum,  and  t  =  0.5,  from 
data  furnished  by  Heath  (9) . 

For  determining  the  thermal  flux  from  the  disintegration  rate,  the 
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expression  is; 


/  Ra«W>F^ 

m  •  Nq»  » (1  ”  e”^^) » (5a 

W=  atomic  mass  of  gold  (197) 

F^  =  ratio  of  thermal  to  total  flux  (0 . 815) 

m  =  mass  of  sample 

JSTq  =  Avogadro's  number  (6.02  x  10^^) 

Ffd  =  flux  depression  correction  (0.99) 

Q-'Xt  =  correction  for  decay  of  sample 
1  -  e”^'^  =  activation  factor 

6'a  =  effective  cross  section  of  sample  (121.5  barns) 

was  calculated  from  the  cadmium  ratio  of  5.36  determined  by 
Kelly  and  Clements. 

Ffij  was  calculated  for  a  0.0005"  foil  by  the  method  of  Ritchie  and 
Eldridge  (12) . 

(J'a  was  calculated  by  applying  the  spectral  hardening  effect  found 
by  Cooke  (3)  to  the  procedure  developed  by  Westcott  (14)’  for  deriving 
effective  from  thermal  cross-sections.  The  basic  expression  is; 

6^  a  =  6o  (g  +  rs) 

5'a  =  effective  cross-section  (121.5  barns) 

<5o  =  thermal  cross-section  (98.8  barns) 
g  F  non  -  l/v  factor  (1.0053  at  20® C) 

r  =  an  expression  relating  thermal  to  total  flux  (0,013  for  this 
reactor) 
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s  =  correction  for  resonance  capture  (17 » 3  at  20  C) 

From  the  data  in  Westcott’s  paper  and  the  cadmium  ratio^  r  was 
calculated^  and  the  values  found  used  to  determine  the  effective  capture 
cross-section,  which  is  higher  than  the  thermal  because  of  the  large 
resonance  correction  for  gold. 

The  results  derived  from  the  experimental  data  are  presented  in  Table 
I,  and  a  curve  of  flux  versus  nominal  power  is  given  in  Figure  3. 


Power 

(watts) 

TABbE  I 

Total  flux 
(onVcmVsec) 

Thermal  flux 
(onVcmVsec) 

Ool 

7.13  X  10® 

(5.81  +  .12)  X  10® 

0«1 

7o29  X  10® 

(5o95  ±  .12)  X  10® 

1 

6.73  x  l07 

(5.48  ±  .11)  X  10^ 

1 

7.02  X  10^ 

(5.72  +  .11)  xl07 

10 

6.44  X  10® 

(5.25  +  .10)  xlO® 

10 

6o41 X  10® 

(5.22  +  .10)  X  10® 

100 

6.45x10® 

(5.26+  .10)  X  10® 

100 

6o50  X  10® 

(5.30  +  .11)  X  10® 

750 

4.61  X  10^° 

(3.75  +  .08)  X  10^® 

750 

4.70x10^° 

(3.83  +  .08)  X  10^® 

8.  Discussion  of  results. 

The  values  obtained  for  thermal  flux  are,  in  general,  higher  than 

those  which  h^ve  been  previously  reported.  Kelly  and  Clements'  figure 

6 

for  Oolw  was  5.31  x  10  ,  while  the  value  given  by  Copeland  and 
Reasonover  was  3.66  x  10®,  and  Swanson  (13)  reported  values  of  5.09  x  10^ 
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7 

4,84  X  10  ,  and  4.32  x  10  for  0,1,  1,  and  10  watts,  respectively. 

His  results  are  especially  interestincf  ^  since  they  were  obtained  from 
good-geometry  beta  counting  of  indium  foils  after  solution  in  HCl.  A 
higher  figure  than  that  of  Copeland  and  Reasonover  is  to  be  expeqted, 
since  they  were  not  aware  of  the  crystal  efficiency  correction  caused 
by  the  air  gap  in  the  crystal  cen;  if  their  figure  is  adjusted  to  the  correct 
efficiency,  and  also  corrected  for  absorption,  it  becomes  5,11  x  10  .  It 
should  also  be  noted  that  between  1961  and  1963  the  physical  location  of 
the  reactor  instrumentation  was  changed.  Although  the  two  control 
channels  affected  were  calibrated  against  each  other  before  and  after  they 
were  moved,  there  is  a  possibility  that  a  particular  indicated  power  level 
is  now  associated  with  a  different  flux. 

The  precision  of  tjie  reported  results  is  essentially  that  of  the 
areas  under  the  photopeaks.  The  average  standard  deviation  of  the  means 
of  all  sets  of  areas  taken  was  2%;  it  is  not  considered  justifiable  to 
adjust  this  figure  to  reflect  the  exapt  deviation  of  a  particular  set  of 
counts,  since  the  deviations  of  all  runs  fell  between  1%  and  3%.  The 
deviations  listed  in  Table  I  are  2%  of  the  reported  flux  values, 

9.  Summary 

Thin  gold  foils  were  activated  in  the  AGN-201  reactor,  and  their 
activation  measured  by  gamma-ray  scintillation  spectrometry.  Areas  of 
photopeaks  were  individually  computed  by  the  GDC  1604  computer  and 
averaged  for  each  run,  in  order  to  minimize  the  effects  of  gain  shift  in 
the  spectrometer.  From  the  measured  photopeak  areas  the  total  and 
thermal  neutron  fluxes  were  computed  for  power  levels  of  0,1,  1,  10,  100, 
and  750  watts.  Flux  values  are  given  in  Table  I.  They  are  somewhat 
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higher  than  the  values  obtained  by  previous  investigators .  Part 
of  the  difference  arises  from  an  incorrect  value  for  crystal  efficiency 
used  in  the  earlier  work . 

I  should  like  to  acknowledge  my  gratitude  and  obligation  to 
Professor  WilUem  W„  Hawes  ^  of  the  Department  of  Metallurgy  and 
Chemistry f  for  guidance,  encouragement,  and  support  throughout  the 
progress  of  this  investigation.  Appreciation  is  also  due  to  Harold  L, 
McFarland,  who  kept  the  reactor  "on  course"  for  many  hours;  tp  Patricia 
C.  Johnson,  who  prepared  the  computer  program  that  is  the  heart  of  the 
calculations,  and  to  Mr.  R.  L.  Heath,  of  Phillips  Petroleum  Co.,  for 
crystal  efficiency  and  other  data . 
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APPENDIX  I 


lo  G-M  coincidence  correction, 

Fpr  corrections  which  are  not  too  large  (less  than  10%  is  the  usual 
criterion) ,  the  relation  of  true  sample  activity  to  the  observed 
activity,  is  given  by: 


1  +  ro(. 


where  X  is  the  dead  time  of  the  Geiger-Muller  tube  „  If  the  counting 
interval  is  chosen  as  unit  time,  the  true  and  observed  total  counts,  n  and 
n  ,  can  be  substituted  for  the  activities ,  In  order  to  determine  X ,  three 
sample^  of  iodine  were  irradiated  until  the  1-128  activity  produced  about  400 
observed  counts  per  second.  Each  sample  was  counted  for  one  minute  at 
six^minute  intervals,  until  the  activity  had  decayed  to  a  level  well  below 
that  for  which  a  correction  would  be  expected,  but  still  high  enough  to 
give  sufficient  statistical  precison.  These  lower  counts  were  thep 
corrected  for  decay  back  to  earlier  times,  giving  a  series  of  computed 
"true-'  counts,  n^. 

Substituting  and  rearranging  in  the  equation  above  gives: 

nt  n't 

The  observed  count  for  each  of  several  times  v^as  used  with  the 
computed  nt  for  each  sample.  Statistical  variation  in  the  counts,  which 
is  multiplied  when  the  decay  correction  is  applied,  introduces  variations 
in  the  computed  value  of  Ve  but  by  taking  enough  values  a  good  mean  can 
be  obtained,  In  the  present  case,  seven  values  were  chosen  from  each 
of  three  samples.  Values  for  X computed  from  the  first  few  counts  are  not 
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usable,  either  because  the  coincidence  correction  is  greater  than  10% 
or  because  the  resolution  time  of  the  scaler  becomes  important,  but 
with  a  colugin  of  computed  values  it  is  easy  to  see  the  point  a|  |^ich 
they  settle  down  around  a  steady  mean „  A  similar  precaution  applies 
to  using  counts  below  the  middle  range.  If  the  statistical  variation  of 
n^t  is  too  great,  V  will  be  unreliable . 

2  o  Initial  activity  from  integrated  counts . 

In  the  present  work,  the  purpose  cf  counting  the  indium  foils  was  to 
determine  their  activity ,0(0,  at  the  time  of  removal  frpm  the  reactor.  In 
order  to  have  the  best  precision,  it  was  desirable  to  count  the  sample  over 
a  fairly  long  period  of  time  —  10  minutes  was  chosen  —  but  with  its  54- 
minute  half-life,  In^ll6  decays  significantly  in  10  minutes,  and  both  the 
aetivity  and  the  coincidence  correction  change,  at  different  rates,  be¬ 
tween  the  beginning  and  end  of  the  count.  The  following  derivation 
enables  one  to  calculate  ogdirectly  from  an  integrated  count  taken  oyer 
any  period. 

Let 

XH  =  true  activity  at  any  time . 

Of'^=;  observed  activity  at  any  time. 

<^o~  he  sired  activity  at  time  to. 
o<t  =  activity  at  start  of  count, 
tj  =?  time  from  to  to  start  of  count. 

T  =  length  of  count. 

'I:' =  G-M  dead  time , 
n  =  true  integrated  count  in  time  T. 
n^=  observed  integrated  count  in  time  T. 


Z2 


By  definition  i 


.T 


And 


n  =  /Cy  dt 

n'  =  /«'  dt  =  f  °*  * 

/  1  +  M 

'0 


dt 


n - dt 


0(idt 


0(i 


1  +  roi  I 


Xo^i  +  e>‘t 


- -  dt  =  0(1 


,r 

dt 


l  +  W 


0  Voll  +  e>“ 


r 


n  = 


di 


OL^T  0(i>vT 


1  In  (0(iT+ 


T 


/  T 

n  = 


r 


'  1 

In  (o(j:  7r  +  + 


>\r 


n/=  T  1  O^i^  +  1 

^  0({t+  e 


fY  •  O"  "t  1  /  T 

Then  In-  ^  _  =  (n^  - -T:r-i>lY=  n^T-Xl 


C/iTr+  \  r 


In  (OfiT  t  1) 
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In  ej^ponential  form, 

-  en>t->,T  e 

Solving  for  i , 

-  1 

^  ^  _ _ ^ _ 

^  r(l  -  e->‘'I) 


From  which 

(^n  1)  g>\t, 

(1  _ 
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APPENDIX  II 


EXPERIMENTAL  DATA 

Experimental  points  us^d  in  determining  photopeak  ar^as  are  given 

LeloWo  Each  pet  of  fqur  points  represents  one  complete  count  over  the 

photopeak;  all  counts  were  one  minute  long.  Valuqs  givep  are  corrected 

back  to  the  tifne  of  removal  from  the  reactor. 

0.1  watt  Run  1 

Mass;  30.0  mg 

Irracjiation  time  477.7  min. 


.Channel 

Count 

1.  19 

3020 

20 

2778 

31 

1769 

22 

1038 

2,  19 

2993 

20 

2917 

21 

1729 

22 

1136 

3.  19 

3496 

20 

2889 

21 

1625 

22 

752 

4.19 

3103 

20 

2759 

21 

1510 

22 

936 

Channel 

Count 

5o 

19 

3093 

20 

2710 

21 

1676 

22 

832 

6. 

19 

3|19 

20 

2732 

21 

1540 

22 

855 

7. 

19 

3107 

20 

2707 

21 

1801 

22 

875 

8 . 

19 

3050 

20 

2790 

21 

1776 

22 

951 

0. 1  watt  Run  2 

Mags:  31.5  mg 

Irradiation  time  534  min. 


19 

5161 

3o 

19 

5337 

20 

5408 

20 

5315 

21 

4191 

21 

3628 

22 

2617  • 

22 

2318 

2 .  19 

5030 

4  o 

19 

5065 

20 

5426 

20 

5343 

21 

4069 

21 

3670 

22 

2443 

22 

2407 

25 


0.1  watt  Run  2  (continued) 


Channel 

Count 

Channel 

5.19 

5016 

6o 

19 

20 

5280 

20 

'21 

4088 

21 

22 

2545 

22 

1  watt 

Run  1 

Ma^s: 

31,0  mg 

Irradiation  time 

60  min. 

1.  23 

5833 

7. 

22 

24 

5340 

23 

25 

3798 

24 

26 

1686 

25 

2.  23 

5861 

al 

22 

24 

5033 

23 

25 

3379 

24 

26 

1711 

25 

3.  23 

5766 

9. 

22 

24 

5187 

23 

25 

3493 

24 

26 

1831 

25 

4.  23 

5965 

10  0 

22 

24 

5020 

23 

25 

3593 

24 

26 

1782 

25 

5 .  23 

6016 

11. 

22 

24 

4623 

23 

25 

3498 

24 

26 

1577 

25 

6.  23 

5855 

24 

4961 

25 

3027 

26 

1368 

Count 

5071 

5204 

4181 

2472 


5276 

5873 

4604 

3043 

5844 

5631 

3957 

2255 

5633 

5647 

4092 

2500 

5725 

5499 

4588 

2997 

5736 

5678 

4514 

2300 


26 


1  watt  Run  2 

Mass;  29 » 5  mg 

Irradiation  time  60  mino 


Channel 

Count 

1,  22 

5204 

5  o 

23 

5359 

24 

4098 

25 

2620 

2.  22 

5349 

6o 

23 

5498 

24 

4347 

25 

2493 

3.  22 

5420 

7, 

23 

5441 

24 

4081 

25 

2499 

4,  22 

5479 

8. 

23 

5522 

24 

4348 

25 

2550 

Chann 

22 

23 

24 

25 


22 

23 

24 

25 


22 

2^ 

24 

25 


22 

23 

24 

25 


10  watt  Run  1 

Mass:  28o5  mg 

Irradiation  time  30  min. 


1.  21 

26739 

6 , 

21 

22 

23023 

22 

23 

14933 

23 

24 

6933 

24 

2.  21 

26929 

7, 

21 

22 

22080 

22 

23 

15087 

23 

24 

7952 

24 

3.  21 

26433 

8. 

21 

22 

23422 

22 

23 

15435 

23 

24 

8062 

24 

4,  21 

26611 

9, 

21 

22 

24276 

22 

23 

16989 

23 

24 

9220 

24 

5.  21 

25863 

22 

24564 

23 

17573 

24 

9287 

Count 

5370 

5295 

3902 

^397 

5444 

5419 

4333 

2509 

5511 

5478 

3222 

2172 

5590 

5374 

3744 

2276 


25665 

24161 

16943 

9^31 

25536 

24034 

16711 

9491 

25144 

23927 

16624 

8849 

26477 

23615 

16310 

6373 
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10  watts  Run  2 

Mass:  30o0  mg 

Irradiation  time  3  0  min . 


Channel 

Count 

Channel 

Count 

1,21 

27346 

5. 

21 

27427 

22 

24734 

22 

24370 

23 

mi\ 

23 

16329 

24 

7871 

24 

8544 

2.  21 

27ii0 

6. 

21 

27886 

22 

24814 

22 

23  8?  7 

23 

17601 

23 

15840 

24 

9104 

24 

7741 

3.  21 

27999 

7o 

21 

27982 

22 

24786 

22 

24177 

23 

16980 

23 

16002 

24 

8022 

24 

7915 

4,21 

27271 

8. 

21 

27862 

22 

24834 

22 

23279 

23 

16740 

23 

15290 

24 

9071 

24 

7799 

100  watts  Run  1 

Mass:  29.0  mg 

Irradiation  time  10  min. 


i.  2P 

90328 

6. 

20 

90152 

21 

79321 

21 

78586 

22 

54308 

22 

53088 

23 

25996 

23 

25587 

2.  20 

92303 

7. 

20 

85047 

21 

82679 

21 

85918 

22 

50305 

22 

65734 

23 

24699 

23 

38091 

3.20 

89959 

8. 

20 

86027 

21 

80298 

21 

86495 

22 

54834 

22 

62329 

23 

27463 

23 

33931 

4.  20 

91006 

9 . 

20 

89215 

21 

74540 

21 

83291 

22 

45774 

22 

58807 

23 

21665 

23 

32023 

5.  20 

89693 

10. 

20 

90518 

21 

79115 

21 

79197 

22 

52516 

22 

51580 

23 

27593 

23 

26720 
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100  watts  Run  2 

Mass:  31.0  pig 

Irradiation  time  10  min. 


Channel 

Count 

Channel 

Count 

1. 

97092 

6. 

20 

97500 

21 

83991 

21 

86819 

22 

57262 

22 

57588 

23 

29511 

23 

28552 

2.  20 

95656 

7o 

20 

95811 

21 

82661 

21 

85773 

22 

53035 

22 

59252 

23 

27731 

23 

27300 

3.  20 

96532 

20 

97178 

21 

81315 

21 

84094 

22 

53139 

.  22 

57617 

23 

24804 

23 

30875 

4c  20 

95844 

9, 

20 

95279 

21 

84405 

21 

89158 

22 

55160 

22 

60150 

23 

28091 

23 

30561 

2Q 

96219 

'  -m 

80541 

.  tz 

54445 

> 

23 

25940 

750  watts 

Run  1 

Mass: 

30,0  mg 

Irradiation  time 

5  min. 

1.  23 

361,355 

5 . 

23 

353,481 

24 

316,663 

24 

301,519 

25 

206,301 

25 

186,786 

26 

52152 

26 

80748 

2.  23 

360,305 

6c 

23 

360,133 

24 

503,361 

24 

300,968 

25 

202,571 

25 

197,930 

26 

91250 

26 

91893 

3.  23 

358,049 

7. 

23 

363,002 

24  , 

306,995 

24 

316,068 

25 

20p,594 

25 

210,761 

26 

88118 

26 

93839 

4.  23 

350,035 

24 

293,764 

25 

186,034 

26 

77688 
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750  watts 

Run  2 

Mass: 

32 . 0  mg 

Irradiation  time 

5  min. 

Channel 

Count 

Channel 

Count 

1.  23 

387,878 

1-  ■ 

23 

381,343 

24 

342, 06p 

34 

359,159 

25 

222,754 

25 

221,929 

26 

105,531 

26 

101,741 

2.  23 

383,478 

5o 

23 

392,784 

24 

340,082 

24 

341,144 

2S 

212,931 

25 

216,476 

26 

99094 

26 

102,292 

3o  23 

387,472 

24 

335,694 

25 

208,157 

26 

96726 

APPENDIX  III 
PHOTOPEAK  PAR;\METERS 


Abscissa  of 

Ordinate 

Foil 

Count 

maximum 

of  maximum 

Area 

0.1  Wj  Run  1 

1. 

18.93 

4969 

25753 

2r 

18o91 

4992 

26860 

3., 

18o76 

5772 

26513 

4.. 

18o43 

5340 

29662 

5. 

18„97 

5052 

23569 

6. 

18o67 

5223 

26492 

7. 

19o05 

5041 

23438 

8. 

19,06 

4997 

23926 

Average 

25777  +  745 

0.1  Wf  Run  2 

1. 

19,66 

5478 

26021 

2, 

19 ,60 

5358 

26082 

3. 

19,21 

5448 

28919 

4. 

19,42 

5273 

26890 

5. 

19,62 

5301 

26104 

6 . 

19,62 

5335 

25  812 

Average 

26638  +  480 

1 W/  Run  1 

1. 

23.36 

6913 

24872 

2 . 

?2.94 

5855 

28620 

3. 

23,08 

5781 

27824 

4. 

22,97 

5915 

29148 

6. 

22,81 

5924 

29585  , 

6 . 

22,96 

5863 

26170 

7. 

22,78 

5765 

28144 

8. 

22,29 

5937 

28885 

9. 

22,40 

5787 

28942 

10, 

22,39 

5782 

33146 

11. 

22,60 

6017 

26367 

Average 

28337  +  592 

1  Ws  Run  2 

1. 

22,53 

5413 

27661 

2. 

22,63 

5653 

26287 

3. 

22,45 

5581 

28033 

4. 

22,56 

5705 

27530 

5« 

22.36 

5480 

2  8094 

6. 

22.55 

5647 

27295 

7. 

21,95 

5671 

30622 

8, 

22.19 

5660 

^9339 

Average 

28107  +  469 
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Foil 

10  w;  Run  1 


10  Wf  Run  ? 


100  w,  Run  1 


100  w.  Run  2 


AjDScissa  of 

Ordinate 

Count 

maximum 

of  maximum 

Area 

lo 

21o04 

26692 

120.798 

2o 

20o85 

27201 

145.846 

3. 

21.01 

26504 

128  .677 

4, 

2I0I4 

26691 

131,213 

5o 

21o34 

26271 

121., 741 

6„ 

21o26 

25952 

122,990 

7. 

2I0I9 

25769 

127,843 

8. 

21.30 

25532 

118,468 

9o 

21.35 

26788 

105.950 

Average 

124.836  +  3599 

Ip 

21o29 

27545 

118.689 

2o 

21o23 

27236 

128,184 

3, 

2I0I7 

27993 

125.965 

4» 

21o09 

27394 

134.143 

5. 

2I0O3 

27471 

133,548 

6  0 

20o97 

28103 

140,952 

7o 

20o99 

27954 

132.878 

8p 

20o78 

28016 

141.493 

Average 

130,914  +  2378 

lo 

20.13 

90149 

412,613 

2o 

20.05 

92989 

420.643 

3. 

20.13 

90021 

420,746 

4p 

19.76 

91841 

438,778 

5p 

19o98 

89884 

441,322 

6p 

2O0O8 

89932 

416,033 

7o 

20o54 

88519 

420,892 

8p 

20o49 

89438 

403,893 

9« 

20o23 

89946 

433.832 

10  „ 

19o95 

90793 

442.734 

Average 

425,149  +  4192 

lo 

19o97 

96940 

478,690 

2o 

19  086 

96248 

478,386 

3o 

l^o96 

96233 

446,182 

4o 

2O0OI 

96038 

458,672 

5o 

19o94 

95733 

455,721 

"6o 

20.10 

97662 

452,845 

7o 

20o24 

98128 

421,151 

8. 

19o90 

97235 

498.277 

9. 

20o26 

96441 

436.508 

Average 


458.492  +  7866 
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750  w.  Run  1 


75 0  w.  Run  2 


Abscissa  of 

Ordinate 

Count 

maximum 

of  maximum 

Area 

1. 

23ol5 

36i;455 

1.367.600 

2o 

23o03 

358.013 

1.619.320 

3„ 

23.10 

356.814 

1.557.260 

4. 

23o06 

348.365 

1.489.360 

5o 

23o05 

353.070 

1.521.690 

60 

22o95 

358.724 

1.669.420 

7o 

23.15 

362.365 

1.579.730 

Average 

1.572.050 

lo 

23.10 

388.279 

1.747.580 

2o 

23.10 

384.938 

1.694.390 

3o 

23.01 

388,204 

1,738,740 

4« 

23.28 

388.743 

1.613.890 

5o 

23.03 

393.131 

1.784.240 

Average  1,715^770 


22570 


29210 


33 


